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Stable oscillation in spin torque oscillator excited by a small in-plane magnetic field 
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Theoretical conditions to excite self-oscillation in a spin torque oscillator consisting of a perpendic¬ 
ularly magnetized free layer and an in-plane magnetized pinned layer are investigated by analytically 
solving the Landau-Lifshitz-Gilbert equation. The analytical relation between the current and oscil¬ 
lation frequency is derived. It is found that a large amplitude oscillation can be excited by applying 
a small field pointing to the direction anti-parallel to the magnetization of the pinned layer. The 
validity of the analytical results is confirmed by comparing with numerical simulation, showing good 
agreement especially in a low current region. 


I. INTRODUCTION 

High-frequency devices are in demand for practical ap¬ 
plications such as communication devices, magnetic sen¬ 
sors, and measurement systems. Spin torque oscillator 
(STO) [Ml is a promising candidate for such applica¬ 
tions because high-frequency signal on the order of Gi¬ 
gahertz can be achieved through the giant magnetore¬ 
sistance (GMR) or tunnel magnetoresistance (TMR) ef¬ 
fect. A small size on the order of nano-meter, a wide fre¬ 
quency tunability over than one Gigahertz, and a narrow 
linewidth less than one Megahertz are also preferable fea¬ 
tures for device applications. In particular, an STO con¬ 
sisting of a perpendicularly magnetized free layer (Tsl - fri 
and an in-plane magnetized pinned layer has attracted 
much attention recently [llj | because the large relative 
angle of the magnetizations of the free and pinned layers 
results in a large emission power on the order of 1 /iW 
[HI- Therefore, this type of STO is considered to be the 
model structure for practical applications. 

The self-oscillation of the magnetization in an STO 
is excited when several conditions, such as the direction 
and magnitude of an applied field, are satisfied [igI - I^ . 
While general conditions to excite the self-oscillation had 
been formulated [H, IH, |2^ , the possibility to excite the 
self-oscillation has to be investigated for each problems 
by solving the Landau-Lifshitz-Gilbert (LLG) equation. 
The number of theoretical works studying the conditions 
of the self-oscillation is few because of the complexity 
arising from the nonlinearity of the LLG equation. For 
example, the oscillation properties of STO with the per¬ 
pendicularly magnetized free layer and the in-plane mag¬ 
netized pinned layer were theoretically studied only in the 
presence of a perpendicular field [2lj, in which the axial 
symmetry of the system reduces the difficulty to solve 
the LLG equation. On the other hand, it is still unclear 
whether the self-oscillation can be excited when the field 
direction is changed. It is practically important to inves¬ 
tigate the possibility to excite self-oscillation in the pres¬ 
ence of a field in a different direction due to the following 
reasons. For example in the sensor applications, an exter¬ 


nal magnetic field to be detected points to an arbitrary 
direction. Thus, it is insufficient to study the oscillation 
properties in the presence of a perpendicular field only. 
Another reason relates to the low power consumption. 
Our previous work [2l| showed that the zero-field oscil¬ 
lation is impossible in this type of STO, and a relatively 
high field is necessary to stabilize the self-oscillation by 
a perpendicular field. Discovering a method to excite a 
stable oscillation by a small field is therefore attractive 
from the viewpoint of reducing the power consumption 
of the STO. The change of the field direction might open 
a new possibility to this problem. 


In this paper, we investigate the oscillation properties 
of the STO in the presence of an in-plane magnetic field 
by solving the LLG equation analytically. We find that 
an applied field smaller than Eq. (1161) pointing in the di¬ 
rection anti-parallel to the magnetization of the pinned 
layer enables to excite a stable self-oscillation. This is 
a preferable feature for practical applications, compared 
with the oscillation by a perpendicular field in which rel¬ 
atively large field is required to stabilize the oscillation 
[2lj. The relation between the current density and the 
oscillation frequency is derived. The comparison with nu¬ 
merical simulation guarantees the validity of the theory 
in the low current region. 


The paper is organized as follows. In Sec. m we 
derived the analytical formulas of the work done by 
spin torque and the dissipation due to damping, which 
are necessary to discuss the possibility to excite self¬ 
oscillation, by solving the LLG equation. Using these 
results, we discuss the theoretical conditions to excite 
the self-oscillation by an in-plane magnetic field in Sec. 
uni These analytical discussions are compared with the 
numerical simulation in Sec. m Section El is devoted to 
the discussion. Section lyTl provides our conclusion. 
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FIG. 1: Schematic views of STO under consideration. An ap¬ 
plied magnetic field points to the positive ^-direction, whereas 
the magnetization of the pinned layer is parallel (p = -l-ej;) 
or anti-parallel (p = —e^;) to the applied field. 

II. ANALYTICAL SOLUTION OF LLG 
EQUATION 


The system we consider is schematically shown in Fig. 
[TJ The unit vectors pointing in the magnetization di¬ 
rections of the free and pinned layers are denoted as m 
and p. We assume that the applied magnetic field points 
to the positive x-direction, whereas p is parallel or anti¬ 
parallel to the field. We use the macrospin model in this 
paper because its accuracy was confirmed by recent ex¬ 
periment EH- The LLG equation is given by 

dm n TT , \ dm , , 

—— = — 7 m X H — 'yHsm x (p x m) -|- am x ——, (1) 

(LL (J/L 

where 7 and a are the gyromagiretic ratio and the 
Gilbert damping constant, respectively. In the follow¬ 
ing, we neglect higher order terms of a by assuming that 
a <C 1, which was confirmed experimentally The 

magnetic field H consists of an effective perpendicular 
anisotropy and the in-plane magnetic field iJt as 

H = Hte^ + (2) 

The effective anisotropy is defined as = Hk — 4,71 M, 
where Hk is a crystalline anisotropy while 47 rM is a 
shape anisotropy. The magnitude of the applied field 
is Ht{> 0). The magnetic field is defined as a gradient 
of energy density i? as H = —dE/d{Mm), where 


E = —M Hiirix — 


MH\ 


eff 

K ml. 


( 3 ) 


A spin torque strength is given by Bin 


H. = 


hr]j 


2e(l + Am • p)Md’ 


( 4 ) 


where d is the thickness of the free layer. The spin po¬ 
larization of the pinned layer is 77 . The parameter A de¬ 
termines the dependence of the spin torque strength on 
the relative ang le of the magnetizations, and is usually 
positive |30I432|| . 


The self-oscillation can be excited when the energy 
supplied by the spin torque balances the dissipation due 
to damping, and thereby, the magnetization processes on 
a constant energy curve of E. This condition can be ex¬ 
pressed as -I- Wa = 0, where Ws{E) = 7 M j> dtids[p • 
H — (m • p)(m • H)] and Wa{E) = —a^M j> dt\H^ — (m • 
H)^] are the work done by spin torque and the dissipa¬ 
tion due to damping during a precession on the constant 
energy curve of E, respectively. This averaging technique 
over a constant energy curve was developed by several au¬ 
thors to investigate the magnetization switching rate in 
the thermally activated region and the oscillation proper¬ 
ties in a ferromagnetic multilayer (T^. . In general, 

the magnetization dynamics in the macrospin model is 
described by two variables. The zenith and azimuth an¬ 
gles, or the three components in Cartesian coordinate 
with the constrain |m| = 1 , are often used as the vari¬ 
ables. Instead, energy E and an angle T/f can be chosen 
as the two variables, where ip characterizes the angle of 
the magnetization on a constant energy curve of E [^ . 
The averaging technique over the constant energy curve 
integrates any quantity over 0 < ip < 27r, and reduces 
the number of the variable to one, E. The reduction of 
the number of the variables makes it relatively easy to 
solve the LLG equation. However, there is still a techni¬ 
cal difficulty to analytically investigate the possibility to 
excite self-oscillation arising from the calculations of Ws 
and Wa- To calculate these functions, the solution of the 
precession trajectory of the magnetization on a constant 
energy curve, which is described by dm/dt = — 7 m x H, 
should be derived. In the case of a perpendicular field 
previously studied [ 2 l[, such a trajectory is described by 
trigonometric functions, and the calculations of and 
Wa become easy. On the other hand, in the case of an 
in-plane field, the precession trajectory on the constant 
energy curve is described by elliptic functions. In this 
case, several mathematical formulas are necessary to de¬ 
rive the explicit forms of and Wa- The details of the 
derivation of and Wa in the present system are sum¬ 
marized in Appendix. Here, we show that the explicit 
forms of and are given by 


1f^.{E) 


hrij I ri-rs 


( 5 ) 


'To.i.E) = - 


4aM 


n - rs 


3 Y Hf/(2Ht) 
The functions A4 and Ha are given by 
(1 ± Ari) 


H, 


A4 = T 


1 . 


K(A:)±-E(fc) 
A^fri-ra) A 


I -b A^ ± 2Ari 
A2(ri - r3)(I ± Ara) 


n(fc, n) 


( 6 ) 


( 7 ) 


Ha =Ht 


1 — rl 

ri - rs 


K(k) 


5E 

iU 


3h; 


eff 


HfJ 


E(k), 

( 8 ) 
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where r£ {£ = 1,2,3) are given by Eqs. (IA6|) . (IA7I) . and 
(IMl) . The double sign in Eq. © means the upper when 
p is parallel to the applied field (p = + 62 ,) and the lower 
when p is anti-parallel to the applied field (p = —e^;). 
The first, second, and third kind of complete elliptic inte¬ 
grals are denoted as K(fc), E(fc) and Il{k,n), respectively, 
where 


k = 




r2 - rs 

ri-ra’ 


(9) 


n = 


TA(r2 - rs) 
1 ± Ara 


( 10 ) 


are the modulus and characteristic of the elliptic inte¬ 
gral, respectively; see Appendix. Equation ([5]) is newly 
derived in this paper, while Eq. ([5]) was obtained in Ref. 

M- 

Equations m and m are functions of the energy 
density. The range of the energy density is < 

E < Esaddie, where E^in = -{MHk/ 2)[1 + {Ht/HK)^] 
and Esaddie = —MHt are the minimum and sad¬ 
dle point energies, respectively. The minimum en¬ 
ergy point and the saddle point locate at mndn = 

and rriaaddie = (1,0,0). 
The precession period on the constant energy curve of 
E 


t{E) 


2K(fc) 


( 11 ) 


relates to the oscillation frequency f{E) via 




FIG. 2: The dependence of j{E) on E for (a) p = -j-ej, 
and (b) p = —e^- The values of the applied held Ht are 
5 (solid, red), 50 (dotted, green), and 65 (dashed, blue) Oe, 
respectively. The energy density in the horizontal axis E is 
normalized as [E - Smin)/(Swaddle - i^min) to make Smin = 0 

and F/saddle ~ 1 • 


f{E) = 


1 

V{Ey 


( 12 ) 


III. ANALYTICAL PREDICTION ON 
POSSIBILITY TO EXCITE SELF-OSCILLATION 


Let us discuss the possibility to excite the self¬ 
oscillation in this system. We denote the current density 
satisfying W^E) + Wa{E) = 0 as j{E), 


m = 


AaeMd Ha 
3Hr] Afs 


(13) 


The current j{E) in the limit of E —^ Emin is the critical 
current destabilizing the magnetization in equilibrium, 
which is obtained from the linearized LLG equation as 


Jc = 


2aeMd _ 
hrjV 




eff 


r 1 

{ HtV] 


) _ 


where V is 

^ r (Et/gf) A[1 - \ 

11 ± X{H,/Hf) 2[1 ± J 


(14) 


(15) 


The double sign in Eq. (ITSl) means the upper for p = 
HGx and lower for p = — e^,. We note here for the later 
discussion that jc for p = —e^, changes its sign when the 
applied field magnitude Ht becomes 


Ht = 


(1 - 


K 5 


(16) 


whereas jc for p = +6^ does not. As discussed below, 
the oscillation properties of this type of STO is charac¬ 
terized solely by this critical current This is partic¬ 
ularly different from other type of STO consisting of an 
in-plane magnetized free layer, where two critical cur¬ 
rents are necessary to distinguish the in-plane and out- 
of-plane oscillations [H, 113, Hs, Ii3| • The number of such 
critical currents, as well as their values, is determined by 
the symmetry of the magnetic anisotropy and angular de¬ 
pendence of the spin torque, and should be investigated 
for each system. 

Two conditions should be satisfied to excite self¬ 
oscillation on a constant energy curve of E [ 2 ^ . The first 
one is that j{E) is finite. More practically, j{E) should 
be experimentally available value, which is at maximum 
on the order of 10^ — 10® A/cw? in typical STOs [3-[l,H“ 
ini. The second is that j{E)/jc > 1. This condition 
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relates to the fact that a current larger than jc should 
be applied to excite any kind of magnetization dynamics 
because the magnetization initially stays at the minimum 
energy state. If this condition is not satisfied, the mag¬ 
netization moves to the constant energy curve including 
I^saddie, and eventually stops its dynamics Let us 
study whether these conditions are satisfied for p = ie^,. 
The values of the parameters are taken from the exper¬ 
iments [ll|, where M = 1448 emu/c.c., Hk = 18.6 kOe, 
d = 2 nm, rj — 0.54, A = 77 ^, 7 = 1.732 x 10^ rad/(Oe-s), 
and a = 0.005. The value of Ht is 61 Oe. The depen¬ 
dences of j(E) on E for p = -ke^, and p = —e^ are shown 
in Figs. [2] (a) and (b), respectively. The values of Ht in 
Fig. [2] are chosen as 5, 50, and 65 Oe to show j{E) for 
both Ht < Ht and Ht > Ht- 

First, we discuss the possibility to excite self-oscillation 
when the magnetization of the pinned layer is parallel to 
the applied field (p = -l-ej,). The current j{E) shown 
in Fig. [5] (a) indicates that j{E)/jc < 1 for i? > 
meaning that the second condition mentioned above is 
not satisfied. Therefore, a stable oscillation cannot be 
excited when p is parallel to the applied field. This ar¬ 
gument will be confirmed in Fig. [3] (a) below. 

Next, we discuss the physical interpretation obtained 
from Fig. [3] (b), where the magnetization of the pinned 
layer is anti-parallel to the applied field (p = —Gx). 
When Ht < Ht{= 61 Oe), the condition j{E)/jc > 1 
is satisfied for the region where E is less than a certain 
value, while \j{E)\ diverges at this point. For example, 
j{E) in Fig. (b) diverges at {E - Llmin)/(£^saddie - 
Emin) ^ 0.94 for Ht = 5 Oe and {E - Llmin)/(L^saddie - 
Emin) — 0.23 for Ht = 50 Oe, respectively. These results 
imply that the self-oscillation can be excited by applying 
an electric current larger than |_ 7 c|, in which the preces¬ 
sion amplitude increases with the increase of the current, 
and finally, saturates to a value corresponding to the en¬ 
ergy at the divergence. The oscillation frequency given 
by Eq. m also saturates to a finite value in the large 
current limit. Below, we show that this analytical pre¬ 
diction works well in the low current region, whereas a 
discrepancy between the analytical and numerical results 
appear in the large current region. This discrepancy is 
due to the breakdown of the assumption used in the aver¬ 
aged LLG equation. The divergence of j (E) at a certain 
E means that the net energy supplied by the spin torque 
during a precession, Ws, becomes zero when the magne¬ 
tization arrives at this constant energy curve. Then, the 
spin torque cannot balance the dissipation due to damp¬ 
ing even when a large current is applied, which means 
\j{E)\ —>• 00 . When Ht > Ht, the condition j{E)/jc > 1 
is not satisfied, and thus the self-oscillation cannot be 
excited. 

To summarize the above discussion, the self-oscillation 
can be excited when the magnetization of the pinned 
layer and the applied magnetic field points to the op¬ 
posite directions, and the field magnitude is smaller than 
Eq. (fT31) . When these conditions are satisfied, the ana¬ 
lytical theory predicts that the oscillation frequency con¬ 


tinuously changes from f{Emin) by increasing the cur¬ 
rent density. Equations (HU and m give the relation 
between the current and oscillation frequency. 

IV. NUMERICAL SIMULATION 

We confirm the validity of the above discussion by com¬ 
paring with the numerical simulation. Figure|3](a) shows 
an example of the time evolution of m when p = +ex 
obtained by numerically solving Eq. (HD. The current 
density, 5.0 x 10® A/cm^, is larger than the critical cur¬ 
rent density, 2.2 x 10® A/cm^, estimated by Eq. (fT4ll . As 
shown, the magnetization finally stops at the film-plane, 
and does not show self-oscillation. We confirmed similar 
dynamics when p = —e^,, Ht > Ht, and |j| > |jc|. These 
are consistent with the above discussions. 

On the other hand, a stable self-oscillation is excited 
when p = —ex and Ht{= 5 Oe) < Ht, as shown in Fig. [3] 

(b) , where the current density is 5.0 x 10® A/cm^ while 
the critical current is 2.6 x 10® A/cm^. The dots in Fig. [3] 

(c) are the peak frequencies of the Fourier transformation 
of rrix, which is measured in experiment [llj |. The ana¬ 
lytical estimation from Eqs. (USD and (USD, shown by the 
solid line in Fig. |3] (c), works well in the relatively low 
current region, guaranteeing the validity of the theory. 
On the other hand, a discrepancy between the analyt¬ 
ical and numerical results is found in the high current 
region due to the following reason. The theory assumes 
that the averaged work done by spin torque and dissi¬ 
pation due to damping during a precession are equal in 
the self-oscillation state. At each point on the constant 
energy curve however, the magnitudes of the spin torque 
and the damping torque are not equal because the spin 
torque and the damping torque have different angular 
dependences. When the spin torque (current) magni¬ 
tude or the damping constant is large, a shift from the 
constant energy curve becomes large, and thus the mag¬ 
netization cannot return to the constant energy curve 
during a precession. Then, the assumption in the aver¬ 
aged LLG equation does not stand, and the analytical 
theory is no more applicable to reproduce the numeri¬ 
cal results. In fact, when the current density is 10 x 10® 
A/cm^, the magnetization does not show self-oscillation. 
In stead, it moves to direction, and stops its dynamics 
because of the large spin torque strength. Therefore, the 
numerically obtained peak frequency in Fig. [3](c) is zero. 

V. DISCUSSION 

The above analytical and numerical results show that 
a stable self-oscillation can be excited when a magnetic 
field smaller than Eq. (fT31) is applied anti-parallel to the 
magnetization of the pinned layer. The oscillation fre¬ 
quency changes continuously by increasing the current 
magnitude. The analytical relation between the current 
and frequency, obtained from Eqs. dUD and m, works 
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FIG. 3: The time evolutions of m at y = 5.0 x 10® A/cm^ with Ht = 5 Oe for (a) p = +ea; and (b) p = — e^;. The inset 
in (b) is an enlarged view, (c) The dependence of the oscillation frequency for p = — e^: on the current density. The dots are 
obtained by the numerical simulation, while the solid line is obtained from Eqs. m and US. It should be noted that these 
equations provide the relation between f{E) and j{E) above jc, which is 2.6 x 10® A/cm^ for Ht = 5 Oe. 


well in the low current region. In a high current region, 
the magnetization eventually arrives at the film-plane, 
and stops its dynamics. In this section, we compare these 
results with the previous works, and discuss the advan¬ 
tages to use the in-plane held. 

Let us compare the present results with our previous 
work (^ . where the oscillation properties in the pres¬ 
ence of a perpendicular held along the z-axis were stud¬ 
ied. It was shown that, when the perpendicular held is 
smaller than He = [3A^/(2 — 3A^)]iL^®, the magnetiza¬ 
tion discontinuously moves from the perpendicular di¬ 
rection to a large cone angle by applying a current larger 
than the critical current. This discontinuous change of 
the cone angle signihcantly reduces the range of oscilla¬ 
tion frequency in the self-oscillation state, as conhrmed 
by experiment [4l[. Thus, a held larger than He should 
be applied to observe the oscillation in a wide range of 
the frequency [ll|, [2lj . On the other hand, in the case 
of the in-plane held studied here, the cone angle of the 
magnetization continuously changes with the increase of 
the current. Therefore, a wide range of the oscillation 
frequency can be used. Also, the self-oscillation can be 
excited by a small held less than Ht in Eq. (USD, which is 
preferable from the view point of a low-power consump¬ 
tion. Therefore, the use of an in-plane magnetic held will 
be an important technical tool for piratical application. 

The system studied above looks similar to an STO with 
the held-like torque (43 - l43| . which was studied in Ref. 
[i^l , because both the held-like torque and the torque due 
to the in-plane held point to m x e^, direction. However, 
there are important differences between these torques. 
The strength of the held-like torque depends on the elec¬ 
tric current, and thus, the spin torque and the held-like 
torque cannot be manipulated independently. On the 
other hand, the magnitude of the in-plane held can be 
controlled independently from the spin torque. The sign 
of the held-like torque is determined by the material pa¬ 
rameters and system size [42 - 143 . 141ll - l43| , while the direc¬ 
tion of the in-plane held is controllable. Regarding these 
points, the in-plane held is more convenient than the 


held-like torque to control self-oscillation, although the 
use of the held-like torque enables to excite self-oscillation 
without an external held. 

Finally, we point out the importance of the parameter 
A. The angular dependence of the spin torque, Eq. 0, 
characterized by A, arises from the spin-dependent trans¬ 
port in a ferromagnetic multilayer. As can be seen from 
Eq. ([7]), the spin torque strength near the anti-parallel 
alignment of the magnetizations (m ~ — p) is larger than 
that near the parallel alignment (m ~ p). This is be¬ 
cause, for example in a GMR system, the amount of the 
spin accumulation, which gives the spin torque on the 
magnetization, is large when the relative angle of the 
magnetizations is large. This parameter had been often 
neglected in previous works on, for example, an in-plane 
magnetized STO [1^ because this angular dependence 
gives small corrections on the oscillation properties such 
as the critical current. However, the parameter A plays 
the dominant role in STO in which the constant energy 
curve is symmetric with respect to an axis perpendicu¬ 
lar to the pinned layer magnetization [^ . as in the case 
of the present STO with a perpendicular field [2l|. Al¬ 
though the constant energy curve in the this paper is 
asymmetric due to the in-plane field, this parameter still 
plays a key role to excite the oscillation. This is because 
the held Ht satishes the following relation; 


lim Ht = 0. (17) 

A^O 


This indicates that the self-oscillation cannot be excited 
by an in-plane held when A = 0. In other words, a hnite 
A, which characterizes the asymmetry of the spin torque, 
is necessary to excite the self-oscillation. We notice that 
this conclusion is consistent with the previous study on 
an impossibility to excite self-oscillation by spin Hall ef¬ 
fect 2^, where the spin Hall system corresponds to the 
limit of A —>■ 0. 
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VI. CONCLUSION 

In summary, it was shown that a stable self-oscillation 
can be excited in an STO with a perpendicularly magne¬ 
tized free layer and an in-plane magnetized pinned layer 
by applying a magnetic field smaller than a critical value, 
Eq. (1161) . pointing in the direction anti-parallel to the 
magnetization in the pinned layer. This is a preferable 
feature compared with the self-oscillation by a perpen¬ 
dicular field, in which relatively large field is necessary 
to stabilize the self-oscillation. The analytical relation 
between the current and the oscillation frequency was 
derived by solving the LLG equation. The comparison 
with numerical simulation shows that the theory works 
well in a low current region. 
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Appendix A: Derivation of Eqs. ([5]) and ([6]) 

The functions Ws{E) and generally defined as 
W^{E) =^M (f dtHs [p • H - (m • p) (m • H)] , (Al) 


^c.{E) 


—a^M 


dt 




(m • H)^ 


(A2) 


are the work done by spin torque and the dissipation 
due to damping during a precession on a constant energy 
curve of E. The integral is over a precession period, 
given by Eq. (HU, on the constant energy curve of E. 
Thus, the precession trajectory of the magnetization on 
the constant energy curve should be substituted into m 
in Eqs. (IA1| and (IA2I) . In the present system, such m is 
given by [24| 

mx{E) = ra -I- (r 2 - r 3 )sn^(it, fc), (A3) 


my{E) = —(r 2 — r 3 )sn(u, fc)cn(u, k) (A4) 


Wz(E) = \Jl-rl-{rl- r|)sn2(u,fc), (A5) 

where u = EltEl^/2y/ri — r^t, and re are given by 


r2{E) 


Ht 

ReS 



Ht 

Hf 


2 


2E 

MHf’ 


(A7) 


rsiE) 


Ht 

Hf 



2E 

MHf’ 


(A 8 ) 


The Jacobi elliptic functions, sn(u, fc), cn(u, fc), and 
dn(M, k) are useful to describe the precession trajectory. 
In the absence of the in-plane field, Ht = 0, the sys¬ 
tem shows the axial symmetry along the z-axis. In this 
case, the modulus k becomes zero, and thus, the elliptic 
functions become trigonometric functions. Then, m in 
Eqs. (TAsI) . (lAI)) . and (IA5I) become = — sin^coswt, 
my = — sin^sinwf, and m^ = cos 0 , where 0 is a con¬ 
stant and oj = "fH^ cos 9. This solution describes circu¬ 
lar motion of the magnetization around the z-axis with 
the initial condition mx(0) = —sinO. 

Substituting Eqs. (IA3I) . (IA4L and (IA5|) into Eqs. (ED 
and (jA2p . we obtain Eqs. m and (ED- The explicit form 
of the integral of Eq. ED for p = 4 - 63 ; is given by 


r f “ ^^ 3 ) - Hf[r 2 +r 3 + A(1 - r|)] 
^ 2 ed Jo { A 2 

Ht{r2-r3)-Hf{rl-rl) ^ 

---sn [u^k) 

A 

_ Ht{l - A^) - gjf [r 2 +r3 + A(1 -f r 2 r 3 )] | 
A^[l -I- Xr3 + A(r2 - r3)sn^{u, fc)] J ’ 

(A9) 


while that of Eq. ([BD is shown in Ref. 2^ . The explicit 
form of Ws for p = —is obtained by replacing A with 
—A and multiplying a minus sign in whole. Then, Eqs. 
([S]) and dp are derived by using the following integral 
formulas 1 ^: 


F[am(u, k), k] — E[am(u, fc), k] 


/ du'sv?{u',k) = ^2 

/■" , 4 . , sn(u, fc)cn(u, fc)dn(M, fc) 

/ du sn [u , k) =- 


3fc2 
2 J- 

+ 3fc4 F[ain(M, fc),fc] 

2(1 + ^ n 

-—E[am(u, fc), fcj. 


(AlO) 


(All) 


/ u /-7 / 1 

— - ^ = ■ 7 n[am(M, k), fc, -a/b], (A12) 

0 -I-asn"'(it', fc) b 


where 


f{v,k)= f =, (A13) 

^0 V 1 — sin 


riiE) = 


E 

Mi7t’ 


(A 6 ) 


E{^,k) 




(A14) 
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Jo (1 — nsin^ T’O'v/l ~ sin^ (p' 


(A15) 


are the first, second, and third kind of incomplete elliptic 
integrals, respectively. The Jacobi amplitude function is 
denoted as am(M, A:). We note that F[am(M, A:), fc] = u, 


E[am(K(A:), fc). A:] = E(fc), and n[am(K(A:), fc). A:, n] = 
n(fc,n), where the first, second, and third kind of com¬ 
plete integrals are defined as K(A:) = F(7r/2, A:), E(fc) = 
E(7r/2,A:), and n(fc,n) = n(7r/2. A:, n), respectively. 
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